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Abstract 

A software platform for global optimisation, 
called PaGMO, has been developed within the Ad- 
vanced Concepts Team (ACT) at the European 
Space Agency, and was recently released as an 
open-source project. 

PaGMO is built to tackle high-dimensional 
global optimisation problems, and it has been suc- 
cessfully used to find solutions to real-life engi- 
neering problems among which the preliminary de- 
sign of interplanetary spacecraft trajectories - both 
chemical (including multiple flybys and deep-space 
maneuvers) and low-thrust (limited, at the mo- 
ment, to single phase trajectories), the inverse de- 
sign of nano-structured radiators and the design of 
non-reactive controllers for planetary rovers. 

Featuring an arsenal of global and local opti- 
misation algorithms (including genetic algorithms, 
differential evolution, simulated annealing, parti- 
cle swarm optimisation, compass search, improved 
harmony search, and various interfaces to libraries 
for local optimisation such as SNOPT, IPOPT, 
GSL and NLopt), PaGMO is at its core a C-l-l- 
library which employs an object-oriented architec- 
ture providing a clean and easily-extensible opti- 
misation framework. Adoption of multi-threaded 
programming ensures the efficient exploitation of 
modern multi-core architectures and allows for 
a straightforward implementation of the island 
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model paradigm, in which multiple populations of 
candidate solutions asynchronously exchange infor- 
mation in order to speed-up and improve the op- 
timisation process. In addition to the C-|— I- in- 
terface, PaGMO's capabilities are exposed to the 
high-level language Python, so that it is possible 
to easily use PaGMO in an interactive session and 
take advantage of the numerous scientific Python 
libraries available. 



1 Introduction 

With the introduction of mass-produced multi-core ar- 
chitectures, personal computers are becoming increas- 
ingly capable of performing parallel computations. Yet, 
the effort to parallelize algorithms is time consuming 
and often not attractive, especially in scientific com- 
puting where software reuse is not as spread a prac- 
tice as in other fields of computing. The open-source 
project PaGMO, (Parallel Global Multiobjective Op- 
timiser), aims at filling this gap for optimisation al- 
gorithms providing, through a generalization of the so 
called island model (i.e. a coarse grained approach to 
parallelization of genetic algorithms) to all types of algo- 
rithms (population based and not), a simple experiment- 
ing platform that allows scientists to easily code algo- 
rithms and problems without having to care at all about 
the underlying parallelization that is provided 'for free' 
by the PaGMO infrastructure . The resulting software 
platform, participating to the Google initiative Summer 
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of Code 2010, is described in this paper together with 
apphcation examples to real life engineering problems of 
interest to aerospace engineers. 



Recent results in global optimisation algorithms ap- 
plied to the design of chemically-propelled interplane- 
tary trajectories have shown how a straightforward ap- 
plication of off-the-shelf optimisation algorithms does 
not suffice to find satisfactory solutions for the most 
complex cases such as the Messenger trajectory or the 
Cassini or the TandEM trajectory [15]. While a wise 
use of the algorithms still provides useful information 
also in these most complex cases, the final optimal so- 
lutions need a substantial amount of engineering knowl- 
edge to be found. In this paper we show how the use of 
PaGMO allows different algorithms to cooperate to the 
solution of the same interplanetary trajectory problem, 
allowing to find solutions also in the most difficult cases 
in a reasonable time. In particular, we demonstrate a 
fully automated search of the solution space based on 
the use of Differential Evolution, Simulated Annealing 
and local search in a cooperative fashion. Information 
is exchanged asynchronously between the solvers oper- 
ating in parallel CPUs via the implementation of a gen- 
eralized migration operator offered by PaGMO. We test 
this search strategy in the case of the Cassini, TandEM 
and Messenger trajectories as defined in the European 
Space Agency Global Trajectory optimisation Problems 
database (GTOP) [17, 44]. We show that the algorithms 
are able to locate interesting regions of the search space 
and in particular to find the possible resonances. In 
the case of TandEM and Cassini, an automatic pruning 
strategy is able to successfully identify the best known 
solutions. In the case of Messenger, the automated 
search locates a large number of possible solution clus- 
ters (due to the possible resonances at Mercury) . A sec- 
ond run of the search focussed on a particular one of 
these clusters holds satisfactory results and is, in partic- 
ular, able to find the same strategy adopted by the actual 
Messenger mission. A last example is then presented, 
where 4406 simpler interplanetary trajectories are op- 
timised (each five times) taking dvantage of PaGMO's 
parallelization capabilities in a reasonably short time to 
locate preliminarly good targets for an asteroid sample 
return mission in the 2020-2050 time frame. 



2 PaGMO 

PaGMO is an optimisation framework developed within 
the Advanced Concepts Team of the European Space 
Agency. Written in C-I--I-, PaGMO aims to provide an 
extensible infrastructure for defining optimisation prob- 
lems (nonlinear, continuous, integer, mixed-integer, box- 
constrained, nonlinear ly constrained, multi-objective op- 
timisation is supported), coupled with a wide arsenal of 
global and local optimisation algorithms - some of them 
coded directly within PaGMO, others called from exter- 
nal libraries through thin wrappers. At the time of this 
writing, PaGMO provides the following optimisation al- 
gorithms: 

• global and local optimisation algorithms coded di- 
rectly within PaGMO, including a simple genetic 
algorithm [12], differential evolution [43], particle 
swarm optimisation [20], adaptive neighbourhood 
simulated annealing [7], improved harmony search 
[25], compass search [21], monotonic basin hopping 
[46], generalised multistart and Monte Carlo search 
[28]; 

• wrapper for SNOPT [11]; 

• wrapper for IPOPT [45]; 

• wrappers for algorithms from the NLopt library 
[18], including Subplex [35] (an extension of the 
classical Nelder-Mead method), COBYLA [32] and 
BOBYQA [33]; 

• wrappers for algorithms from the GSL library [10], 
including the Broyden-Fletcher-Goldfarb-Shanno 
(BFGS) method [4], Fletcher-Reeves and Polak- 
Ribiere nonlinear conjugate gradient methods [39] 
and the classical Nelder-Mead method [30]; 

• wrappers for algorithms from the SciPy library [19] 
(only available in the Python bindings), includ- 
ing fmin (Nelder-Mead), L-BFGS-B [48], sequential 
least-square programming [22] and truncated New- 
ton method [29]. 

PaGMO provides automatic parallelisation of the op- 
timisation process via a coarse-grained approach based 
on the island model [26], in which multiple optimisa- 
tion instances of the same problem are launched at the 
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same time, asynchronously exchanging information and 
improving the overah convergence properties of the op- 
timisation. In PaGMO's implementation of the island 
model, each optimisation instance (i.e., each island) is 
launched in a separate thread of execution, thus au- 
tomatically taking advantage of modern multiprocessor 
machines. The connections between islands are resp- 
resented by a graph topology in which each node cor- 
responds to an island and the edges represent routes 
through which candidate solutions can be communicated 
from one island to the other. The graph topologies can 
be either constructed by manually adding nodes and 
edges, or they can be selected among those already coded 
within PaGMO, including: 

• popular topologies in the context of parallel 
population-based optimisation, such as fully con- 
nected, torus, cartwheel, lattice, hypercube, broad- 
cast, and various types of ring topologies; 

• small- world network topologies, such as the 
Barabasi- Albert [1] and Watts-Strogatz [47] mod- 
els; 

• G{n,p) Erdos-Renyi random graph [9]; 

• custom topologies (such as the wheel rim topology 
described in §3). 

Some of the topologies available in PaGMO are visu- 
alised in Figure 1. Full control over the fine-grained de- 
tails of the migration strategy (e.g., migration frequency 
and rate, selection and replacement policies) is provided. 
A preliminary study of the impact the topology on the 
optimisation process can be found in [36] . The class that 
contains the set of islands collaborating in an optimisa- 
tion process, the topology and the migration policies is 
known in PaGMO as an archipelago. 

PaGMO ships with a number of implemented optimi- 
sation problems readily available for use, such as: 

• classical continuous test functions, such as Rastri- 
gin, Rosenbrock [34], Schwefel, Griewank, Branin, 
Himmelblau, Lennard- Jones potential [23] and 
Levy 5; 

• constrained continuous test functions from [24]; 

• integer programming problems: Golomb ruler [40], 
0-1 knapsack problem [27]; 



• mult i- objective optimisation test problems from [8]; 

• all the chemical interplanetary spacecraft trajectory 
problems from the European Space Agency's GTOP 
database [17, 44]; 

• an interplanetary multiple gravity assist low-thrust 
problem. 

PaGMO's C-f-|- capabilities are exposed to the high- 
level language Python, so that it is possible to instan- 
tiate problems, algorithms, topologies and islands from 
either a script or an interactive Python session. It is 
also possible to define new problems and algorithms 
directly from Python, thus allowing on one hand to 
rapidly prototype and evaluate new ideas, and on the 
other to leverage the rich ecosystem of freely-available 
scientific Python modules (e.g., numerical integrators, 
machine learning libraries, computer algebra systems, 
etc.). Coupled with the matplotlib plotting module and 
the enhanced Python shell IPython, PaGMO's Python 
bindings (which have been called PyGMO) offer a user- 
friendly interactive graphical experience. 

3 Some examples 

As an example of the use of PaGMO to solve engineering 
problems we report here the results of the application of 
the optimisation strategy described in the previous sec- 
tion to four trajectory optimisation selected problems. 
The first three problems are taken from the European 
Space Agency Global Trajectory optimisation (GTOP) 
database [17, 44]. The problems selected are among 
the most difficult proposed in the database and are in- 
cluded in the basic PaGMO distribution. They all are 
box-constrained, continuous, single objective optimisa- 
tion problems, representing a multiple gravity assist in- 
terplanetary trajectory with one deep space maneuver 
allowed in each trajectory leg. The search space in- 
cludes launch windows spanning decades (see the GTOP 
database for the precise definitions of the allowed bounds 
on the launch and fly-by dates). The fourth problem is 
a simpler problem admitting though a large number of 
different instances. We take advantage of PaGMO par- 
allelization to find solutions to 4406 different instances 
of the problem in a reasonable computing time. 
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Experimental setup All the optimisation problems 
were set up in an archipelago of 5-7 islands (depend- 
ing on the number of available cores on the machine at 
the time of the experiment), equipped with a wheel rim 
topology (see Figure 2). The wheel rim topology con- 
sists of a classical bidirectional ring topology with an 
additional island at the center, fully connected to all the 
other islands. We chose to deploy global optimisation 
algorithms (namely, adaptive neighbourhood simulated 
annealing from [7] and differential evolution from [43]) 
on the ring, and a local optimisation algorithm (namely, 
the Subplex algorithm from [35] as implemented in [18]) 
in the center. 

The motivations behind these choices are the follow- 
ing: 

• the ring topology is a proven and popular choice in 



the context of parallel population-based algorithms, 
as shown for instance in [14, 42, 13, 5, 6, 16, 2]; 

• the additional island in the center receives through 
migration the best results of the global optimisation 
algorithms in the ring, refines them through a local 
search, and migrates them back to the ring. Its 
role is hence, on one hand, to improve the results of 
the global search, and on the other to inject back 
diversified candidate solutions into the ring; 

• regarding the choice of the algorithms, both sim- 
ulated annealing and differential evolution have 
proven to be effective for the optimisation of inter- 
planetary spacecraft trajectories (as shown for in- 
stance in [16]), whereas the derivative- free Subplex 
method, a refinement of the classical Nelder-Mead 
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Figure 2. Experimental setup: an archipelago with 
wheel rim topology, global optimisation algorithms on 
the outer ring (simulated annealing and differential evo- 
lution) and a local optimisation algorithm in the inner 
island (Subplex). 



algorithm, is particularly suited for the noisy and 
multi-modal objective functions appearing in these 
optimisation problems. 

In order to give an example of use of PaGMO, we 
reproduce here the Python code necessary to perform 
one optimisation run with the setup described above: 

1 # Import the PyGMO classes 

2 from PyGMO import * 

3 

4 # Instantiate the algorithms 

5 sa = algorit hm . sa_corana ( 1 00 00 , 1 , . 1 ) 

6 de = algorithm . de (500 ,0.8 ,0.9) 

7 local = algorit hm . n lop t _sbplx ( 500 , 1 e —4) 

8 

9 # Instantiate the problem 

10 prob = problem . messenger_full () 
11 

12 # Build the archipelago 

13 a = archipe lago ( topology . rim {) ) 

14 a.push_back(island( prob , local ,1 ,1.0 , migration. 

worst_r_policy () ) ) 

15 a.push_back(island( prob , sa , 1 ) ) 

16 a.push_back(island( prob , de , 2 ) ) 

17 a.push_back(island( prob , sa , 1 ) ) 



18 a.push_back(island(prob,de,20)) 

19 a. push_back(island (prob ,sa ,1) ) 

20 a.push_back(island{prob,de,20)) 

21 

22 # Perform evolution twenty times 

23 a.evolve{20) 

24 a.join() 

Detailed explanation: 

• on line 2, all the PaGMO classes are imported into 
the current namespace; 

• on lines 5-7, the algorithms are instantiated; 

• on line 10, the problem (in this case the full Mes- 
senger problem) is instantiated; 

• on lines 13-20, the archipelago is instantiated: 

— on line 13, an empty archipelago with rim 
topology is created; 

— on line 14, the central island is created and 
inserted in the archipelago with the push.back 
method. The island is constructed from 
the problem prob and the algorithm local, it 
contains one single individual, has a prob- 
ability of accepting the migrating individu- 
als of 100% and replacement policy migration. 
worst_r_policy 0, which will unconditionally re- 
place the worst individual in the island with 
the incoming individuals. This island needs 
a non-default replacement policy because we 
want it to optimise every candidate solution 
coming from the ring, whereas the default be- 
haviour would be to accept migrating indi- 
viduals only if they would improve upon the 
worst individuals present in the population 
(which is the behaviour frequently desired for 
population-based algorithms); 

— on lines 15-20, the ring islands are created and 
inserted into the archipelago. The simulated 
annealing islands operate on populations of a 
single individual, whereas the differential evo- 
lution islands are instantiated with a popula- 
tion of 20 individuals. The default migration 
policies are in these cases appropriate; 

• on line 23, the optimisation process is started by 
calling the evolve () method of the archipelago. The 
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argument passed to the evolve () method, in this case 
20, means that each algorithm on each island is 
called 20 times with the parameters passed in the 
constructors on lines 5-7. E.g., in case of differen- 
tial evolution, it means that the algorithm is run 
for 500 • 20 = 10000 generations, with weight coef- 
ficient equal to 0.8 and crossover probability equal 
to 0.9. Migration is allowed to happen at the end of 
each one of the 20 internal iterations of the evolve () 
method; 

• on line 24, the archipelago is joined, meaning that 
the flow of the program stops until the optimisa- 
tion run started on line 23 has concluded. Since 
PaGMO runs asynchronously each algorithm in a 
separate thread, the evolve () call on line 23 will re- 
turn almost immediately - the optimisation process 
having forked in the background. The join () call 
blocks the program until the optimisation has fin- 
ished. 

Optimisation strategy For the first three problems 
we adopted the following optimisation strategy: 

1. we instantiated an archipelago with rim topology as 
described above and let the optimisation run for a 
fixed amount of time; 

2. at the end of each optimisation run, we recorded the 
best candidate solution produced and then reset the 
archipelago with randomly-chosen decision vectors. 

Step 1. and 2. where repeated multiple times, thus 
producing a collection of optimised candidate solutions. 
The cluster pruning algorithm described in [15] was then 
run on the collection of candidate solutions, returning 
new problem bounds in which top decision vectors are 
contained. The new bounds are then used to launch 
other rounds of multistart optimisations. 

For the sample return problem, which involves the so- 
lution of different instances of a simpler problem, a single 
run of the optimisation algorithms in the archipelago is 
good enough and thus no cluster pruning was used. 

3.1 Results on problem::cassini_2 

This problem represents the interplanetary trajectory of 
the spacecraft Cassini. For a detailed description of this 



global optimisation problem we refer the reader to the 
GTOP database [17, 44]. For the purpose of this paper 
we just mention that the objective function represents 
the sum of all AV, including the launch, where the last 
contribution (Jupiter arrival) is relative velocity with re- 
spect to Jupiter at arrival. For this problem we apply 
the fully automated search described above with three 
pruning cycles. At the end of the process (employing 
7CPUs for a period of roughly 8 hours) the best tra- 
jectory found is visualized in Figure 3a. Details on the 
trajectory as reported in Table la. The trajectory is, 
essentially, the same best result posted in the database 
(22th May 2009) and found by M. Schliieter, J. Fiala, M. 
Gcrdts at the University of Birmingham using the MI- 
DACO solver developed within the project "Non-linear 
mixed-integer-based Optimisation Technique for Space 
Applications" co-funded by ESA Networking Partner- 
ship Initiative, Astrium Limited (Stevenage, UK) and 
the School of Mathematics, University of Birmingham, 
UK [38]. The solution employs two deep space maneu- 
vers during the first two legs as detailed in Table la. 

3.2 Results on problem: :tandem(6, 10) 

TandEM is one of the L-class candidate missions that 
were proposed in response to the European Space 
Agency call for proposals for the 2015-2025 Cosmic- 
Vision programme. Initially, the mission included a Ti- 
tan orbiter and an Enceladus penetrator. The inter- 
planetary part of the trajectory was preliminarly stud- 
ied in 2008 by a Mission Analysis Outer Planets Work- 
ing Group that included different experts from academia 
and space industry. In that preliminary study a baseline 
for the TandEM mission was defined and forms the basis 
of the problem here solved in an automated fashion us- 
ing PaGMO. The baseline considers a launch with Atlas 
501, and an injection orbit at Saturn with e — 0.985, 
rp ~ 80330 km. The mission objective is to maximize 
the final spacecraft mass at arrival and to complete the 
trajectory within ten years. For a detailed description 
of this global optimisation problem we refer the reader 
to the GTOP database [17, 44]. For this problem we 
apply the fully automated search described above with 
three pruning cycles. At the end of the process (em- 
ploying 7CPUs for a period of roughly 6 hours) the best 
trajectory found is visualized in Figure 3b. Details on 
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Figure 3. Visualisation of the best trajectories found for the first three interplanetary trajectory problems: prob- 
lem:. 'cassini-E (a), problem: :tandem(6, 10) (b) and problernI:messenger-full (c). 
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Table 1. Details of the best trajectories found for the first three interplanetary trajectory problems: problem:: cassini-2 
(a), problem::tandem(6,10) (b) and problem: :messenger_full (c). 



the trajectory are as reported in Table lb. The solution 
found improves the previous best found by B. Addis, 
A. Cassioli, M. Locatelli, F. Schoen (from the Global 
Optimisation Laboratory, University of Florence) who 
also have the record on the solutions for all other Tan- 
dEM problem instances (i.e. for different fly-by and time 
constraint). The final trajectory employs one only deep 
space manouvre during the Venus- Venus leg. 

3.3 Results on problem:: messenger _full 

This problem is probably the most complex problem in 
the GTOP database and represents one of the most com- 
plex chemical interplanetary trajectory ever designed 
and flown, that of the Messenger spacecraft. Messen- 



ger, at the time of writing, is on its way to Mercury, 
where (roughly next year, in 2011) will become the first 
spacecarft to ever orbit around the planet. Its path 
in the solar system to reach its final destination in- 
cluded a long planetary tour: Earth-Earth- Venus- Venus- 
Mercury-Mercury-Mercury-Mercury. In the PaGMO 
version of the problem, the Messenger trajectory is tran- 
scribed into a box-constrained global optimisation prob- 
lem. The objective function is the total AV^ accumu- 
lated from the Earth launch to a Mercury Orbit Inser- 
tion (MOI) into an orbit having e = 0.704, r^ = 2640 
km. For a detailed description of this global optimisa- 
tion problem we refer the reader to the GTOP database 
[17, 44]. In this case, after a first run of the algorithm, 
cluster detection shows the presence of many different 
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Figure 4. Manual pruning for a selection of variables from the Messenger problem. The blue dots represent the best 
candidate solutions of each run of the multistart strategy, the vertical green bars denote the new, narrowed bounds of 
the problem and the red diamond marker represents the final solution. The clustering of the best candidate solutions 
in correspondence with the resonant fiybys at Mercury is clearly visible for the variables T2, T4, T5 and Tq. 
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solution clusters related to the different possible reso- 
nances at Mercury, as shown in Figure 4. The best solu- 
tion after the first algorithm run is around 5.15 km/sec. 
By focussing the optimisation into one of the detected 
clusters we obtain a solution at around 2.3 km/sec which 
is lowering the GTOP database record substantially and 
is detailed in Table Ic and visualized in Figure 3c. 

3.4 Results on problem-samplej-eturn 

A single instance of this problem represents an in- 
terplanetary trajectory starting from the Earth and 
performing a rendezvous with a selected asteroid. After 
a minimum waiting time the spacecraft is required to 
come back to the Earth with a maximum hyperbolic 
encounter velocity. One deep-space maneuver per leg 
was allowed, creating a global optimisation problem 
of dimension 12. This type of trajectory can be used 
to perform the preliminary selection of possible final 
asteroids for sample return missions. The same trajec- 
tory model is also relevant to design human missions 
to asteroids. For the purpose of this paper we do not 
enter into the details on the system design and launch 
window choice, instead it is our interest to 'just pick an 
example' and show the possibility of using PaGMO to 
solve in a reasonable time a large number of problem 
instances (e.g. varying the final asteroid). We took all 
the asteroids listed in the JPL NEA database: 

http : //neo . jpl . nasa . gov/ cgi-bin/neo_elem 

having H < 22, which roughly corresponds to as- 
teroids having diameter larger than 200m. For the 
selected 4406 asteroids we optimised the trajectory con- 
sidering a launch in the 2020-2050 time frame, allowing 
a final encounter velocity at the Earth return of 4.5. 
km/sec and minimizing the total AV as evaluated from: 

AV = AVl + AVdsm, + AVr+ 

AVd + AVds„r2 + ^Ve. 

where AVl is the hyperbolic velocity when leaving the 
Earth sphere of influence, AV^smi is the first deep space 
maneuver, AVr is the rendezvous velocity, AVd is the 
relative velocity at asteroid departure, AVdsm2 is the sec- 
ond deep space maneuver and AVe is the final braking 



maneuver to reduce the entry speed to 4.5 km/sec. A 
minimum waiting time on the asteroid of 5 days is also 
considered. 

The computations where performed on an Xserve with 
8 processing units and lasted 8 hours. The results are 
visualized in Figure 5. 

4 Conclusions and future work 

In this paper we have presented PaGMO, a global op- 
timisation software framework for parallel engineering 
optimisation developed within the Advanced Concepts 
Team at the European Space Agency. We have tested 
PaGMO on three hard, realistic interplanetary trajec- 
tory optimisation problems (TandEM, Gassini and Mes- 
senger), showing how PaGMO is able to find automat- 
ically the best known solutions for all three problems. 
With the help of a human-guided final pruning step, 
PaGMO was also able to locate the 'real' trajectory for 
the Messenger probe, consisting of multiple resonant fly- 
bys at Mercury. A fourth benchmark problem, consist- 
ing in the optimisation of simple trajectories to a selec- 
tion of 4406 near-Earth asteriods, was shown with the 
intent of highlighting PaGMO 's parallelisation capabili- 
ties. 

Future work on PaGMO will concentrate on areas such 
as: 

• extension of the computational capabilities via in- 
terfacing to popular massively parallel frameworks, 
such as MPI [41] for scientific clusters and BOING 
[3] for distributed computing; 

• exploration of the possibility of using GPGPU com- 
puting [31] to speed-up the most time-consuming 
parts of the optimisation process; 

• implementing/interfacing additional optimisation 
algorithms; 

• interfacing with machine learning packages (such as 
PyBrain [37]), for easy coding of artificial intelli- 
gence problems. 

Some of these activities will be tackled within the Google 
Summer of Code 2010, in which an international group 
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of University students will be working during the sum- 
mer on PaGMO under the mentorship of the PaGMO 
development team - while being sponsored by Google. 

PaGMO is Free Software, and it is available for 
download from the SourceForge website: 

http: //pagmo . sourcef orge .net 
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